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We review a systematic practical implementation of the multiple scattering formalism due to 
Balian and Duplantier [R. Balian and B. Duplantier, Ann. Phys. (NY) 104 , 300 (1977); 112 , 165 
(1978)] for the calculation of the Casimir interaction between arbitrarily shaped smooth conductors. 

The leading two-point scattering term of the expansion has a simple compact form, amenable to 
exact or accurate numerical evaluation. It is a general expression which improves upon the proximity 
force and pairwise summation approximations. We show that for many geometries it captures the 
bulk of the interaction effect. The inclusion of terms beyond the two-point approximation provides 
an accuracy check and explains screening. As an illustration of the power and versatility of the 
method we re-evaluate sphere-sphere and sphere-plane interactions and compared the results with 
previous findings that employed different methods. We also compute for the first time interaction of a 
hyperboloid (mimicking an atomic force microscope tip) and a plane. We also analyze the anomalous 
situations involving long cylindrical conductors where the two-point scattering approximation fails 
qualitatively. In such cases analytic summation of the entire scattering series is carried out and a 
topological argument is put forward as an explanation of the result. We give the extension of this 
theory to the case of finite temperatures where the two-point scattering approximation result has a 
simple compact form, also amenable to exact or accurate numerical evaluation. 
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I. MOTIVATION 

Casimir interactions are the macroscopic response of 
the physical vacuum to the introduction of material ob¬ 
jects. They were first derived as an attractive force be¬ 
tween perfectly conductive parallel plates induced by the 
zero-point motion of the electromagnetic field [1] . There 
is convincing experimental evidence for the reality of 
these forces [2] and a vast body of literature dedicated to 
various aspects of the phenomenon [3]. 

The goal of this topical review is to provide a practi¬ 
cal guide to the computation of the Casimir interactions 
between arbitrarily shaped smooth conductors. Since 
the introduction of a conductor into a previously empty 
space modifies the electromagnetic spectrum, and thus 
the zero-point energy, the Casimir energy can be heuris- 
tically presented in the form [T] 

^ = ( 1 ) 

a 

where uja and uja are the mode frequencies before and 
after the introduction of the conductor. This idea can 
be made rigorous by including a smooth cutoff at large 
frequency; a large variety of cutoff functions yield the 
same result [1]. The case of perfect conductors which is 
our focus plays a special role because Casimir forces are 
cutoff-independent thus only reflecting the geometry of 
the problem 

When the conductor is a collection of several disjoined 
objects, part of the Casimir energy depends on the rel¬ 
ative location and orientation of the objects; this is the 
interaction energy Smt- Although it can in principle be 
calculated to arbitary accuracy [MS], in practice such 


calculations encounter significant computational difficul¬ 
ties (the parallel plates geometry is the only exception). 
Then it is useful to have approximate methods of com¬ 
putation, especially when there is a systematic means to 
improve the result. 

Casimir [1] calculated the force between parallel 
planes. A simple method to study non-planar geome¬ 
tries, the Derjaguin or proximity force approximation 
(PEA) extends this to nonplanar objects by treating them 
as a superposition of infinitesimally small facing paral¬ 
lel plane segments [S]. It is assumed that only the fac¬ 
ing sides of the interacting objects need be considered. 
The method is restricted to the case that the distance of 
closest approach is significantly smaller than the objects 
themselves. 

A complementary and equally simple method (pro¬ 
posed by Mostepanenko and Sokolov m) is based on 
the interaction between small objects at large separa¬ 
tion given by Casimir and Polder m, and extending this 
to intermediate sizes and distances by regarding macro¬ 
scopic objects as a collection of infinitesimally small vol¬ 
ume elements and performing a pairwise summation of 
elementary two-body interactions. 

Even though in many geometries these approximation 
schemes give reasonably good results, they suffer from 
two fundamental shortcomings: (i) they assume additiv¬ 
ity of Casimir interactions, though this is known not to 
be the case; (ii) in their original form they are uncon¬ 
trolled, thus making it difficult to independently judge 
and improve the results. 

In recent years significant progress has been made in 
computing analytic corrections to the PEA which is nec¬ 
essary for interpretation of precise modern day measure¬ 
ments of the Casimir force. For example, the leading 
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order proximity corrections (relevant to the case that the 
distance of closest approach is much smaller than object 
sizes) for the interaction between a cylindrical conduc¬ 
tor that is parallel to a conductive plane m or for a 
spherical conductor facing a conductive plane m have 
been computed. These findings nicely fit with the idea 
of an approximation scheme in which the Casimir result 
for parallel plates is supplemented by a phenomenologi¬ 
cal expansion in the transverse derivatives of the height 
function [14) . The approach advocated in this review 
complements these studies because it is not limited by 
the proximity approximation or shape and gives usable 
and improvable approximation at all separation distances 
of the objects. 

This review is based on the multiple scattering expan¬ 
sion technique due to Balian and Duplantier (BD) [5]. 
We will present a systematic approach to the calculation 
of Casimir interactions between arbitrarily shaped con¬ 
ductors, which overcomes the shortcomings of previous 
treatments [gdo]. The interaction energy is represented 
as a sum of terms; the leading term is a two-point in¬ 
tegral, with one point on each surface. We show that 
for many geometries and arbitrary separation of the con¬ 
ductors this term (£2 in what follows) captures the bulk 
of the interaction energy £int^ while being as simple and 
compact as previous approximations iniin]. For these 
geometries the Casimir interactions can be viewed to be 
approximately additive. Including the four-point scatter¬ 
ing contribution £4 gives a useful check on the accuracy: 
usually it is small, but in some problem cases it is not. 

However, there are geometries where the result sup¬ 
plied by £2 is qualitatively wrong in the large separation 
distance limit; in such cases the entire multiple scatter¬ 
ing series is required to obtain the right answer for £int', 
we put forward a topological argument which purports 
to explain such anomalous cases. 

Pertinent to the situations when £2 captures the bulk 
of the effect, we give finite-temperature generalization of 
the result, i.e. the expression for the Casimir interaction 
free energy J -2 which, like £ 2 -, has a simple compact form. 

Despite its generality, as computational technology the 
BD method has received very little attention in the lit¬ 
erature even though a prospect of using it to calculate 
Casimir interactions between arbitrarily shaped conduc¬ 
tors has been mentioned [7]. A 2010 review of the sub¬ 
ject even stated that the BD ’’method has not proved 
workable in practice” m- This is possibly due to com¬ 
plexity of multiple surface integrals involved into evalua¬ 
tion of subsequent terms of the multiple scattering expan¬ 
sion. However these integrals are numerically tractable. 
Specifically we will use Monte Carlo integration which 
allowed us to implement the BD method in the form of a 
flexible practical tool amenable for calculation of Casimir 
interaction between arbitrarily shaped conductors (finite 
element methods would undoubtedly be more efficient 
and precise, but less adaptable). We also note that the 
BD method was recently utilized by us to solve a related 
problem of computation of the Casimir self-energy of an 


arbitrarily-shaped solitary smooth conductor m- There 
is no doubt that the BD method can be also used to 
compute 3- and generally Wbody Casimir interactions 
of arbitrarily shaped conductors. 

The rest of this review is organized as follows: 

In Section H, following BD [^, we give general expres¬ 
sion for the Casimir energy £" as a multiple scattering 
series. 

The part of that energy £int that depends upon the 
relative position and orientation of the conductors is dis¬ 
cussed in Section HI where we observe (Section HIA) that 
the leading two-point contribution of the BD expansion 
£2 has a simple compact form of a double surface integral 
over the two surfaces of the conductors. This expression 
is nearly equivalent to a representation of the Casimir 
interaction as a superposition of elementary Casimir- 
Polder m interactions, which has been proposed on phe¬ 
nomenological grounds [lOj . For simple geometries £2 can 
be computed analytically and several examples of suc¬ 
cessful calculations (two small conductors, small conduc¬ 
tor and a plane, two parallel planes, and two spheres) are 
given. However, there is a qualitative failure of the two- 
point approximation for geometries involving long cylin¬ 
ders. The four-point interaction term is given explicitly 
in Section HIB while higher-order contributions are dis¬ 
cussed in Section HIC. What sets the BD method apart 
from majority of other approaches (outlined in Section 
HID) is that it does not rely on separability of Maxwell’s 
equations which is limited to a few high symmetry cases. 
The only other method known to us where separability is 
not essential is that due to Johnson and collaborators [S] 
which we view as complementary to the BD approach. 

Applications of the BD expansion to a series of con¬ 
crete problems are given in Section IV. We begin with 
the parallel plane geometry (Section IVA) followed by the 
plane-arbitrary surface geometry (Section IVB) where we 
observe that the two-point term of the BD expansion de¬ 
rives (up to numerical factor close to unity) the PFA; a 
counterpart of this result for the sphere-arbitrary surface 
geometry is also given. We then re-evaluate the Casimir 
interaction in the sphere-plane (Section IVC) and sphere- 
sphere (Section IVD) geometries and demonstrate that 
the two-point approximation captures about 90% of the 
result found in the literature while the numerically eval¬ 
uated four-point term accounts for most of the discrep¬ 
ancy. As a potentially important practical application, 
in Section IVE we study the geometry of hyperboloid of 
rotation (mimicking a rounded tip of atomic force micro¬ 
scope) facing a plane, and investigate dependences of the 
Casimir interaction on the plane-tip separation and tip’s 
opening angle. As before, the two-point term is com¬ 
puted analytically while the four-point contribution is 
evaluated numerically. The non-additivity of the Casimir 
effect is discussed in the context of an object near a con¬ 
ducting slab of finite thickness. Here the physical effect 
of screening plays an important role which the two-point 
approximation misses. By analyzing the sphere-slab ge¬ 
ometry (Section IVF) beyond the two-point approxima- 
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tion we explicitly demonstrate that the screening mecha¬ 
nism is built into the BD expansion; simultaneously this 
provides an explanation of the success of the PFA for 
closely positioned conductors. In Section IVG we focus 
on anomalous geometries involving long cylinders where 
the two-point approximation fails qualitatively and the 
entire BD series needs to be accounted for. We relate this 
difficulty to change of space topology due to the presence 
of cylindrical conductor, and then, for the plane-parallel 
cylinder geometry, carry out analytic summation of the 
entire multiple scattering series, and confirm the result 
found in the literature. 

We conclude (Section V) by deriving an expression for 
the Casimir interaction free energy in the two-point ap¬ 
proximation, which has a simple and compact form. 
It encompasses the classical (high-temperature limit) and 
is expected to exhibit the same level of accuracy as its 
zero-temperature counterpart £ 2 - 


II. CASIMIR ENERGY 

BD have given an approach to the Casimir problem 
that allows calculation of £ for a conducting surface of ar¬ 
bitrary shape [16]. The starting point of BD’s method is 
the observation that the electromagnetic response (fields 
and induced currents) to a current source perturbation 
(such as an oscillating magnetic dipole) will diverge at 
the resonant frequencies. Viewed as analytic functions 
of the frequency, the corresponding Green functions have 
poles along the real frequency axis; these can be used to 
construct a contour integral representation for £. This 
is the basis of many approaches to the calculation of the 
Casimir energy dZllISj. 

The magnetic field due to a perturbing current source 
on one side of the surface is completely screened so that 
there is no response to the perturbing current on the 
other. The screening is accomplished by a current jsif) 
on the surface. BD observe that the surface current is 
related by Ampere’s law to the tangential components of 
the total field, which is due both to the perturbing cur¬ 
rent and to the surface current elsewhere on the surface. 
This leads to the integral equation 

]{r2)=3a{r2) + J K{r2,ri) ( 2 ) 

where the integral is over the surface and K is the tensor 
that describes how a current j(fi) at fq would induce a 
current j^rq) at fq, ignoring all the other currents, so that 
the electromagnetic fields propagate as if in empty space. 
The integral equation will give divergent response at the 
mode frequencies. BD show how to use this property 
to calculate the change in the density of states caused 
by the introduction of the conductors. From this the 
Casimir energy can be calculated as a contour integral 
which is most conveniently evaluated by integrating along 
the imaginary wave number axis y = uj/ic. Here the 


tensor K is 

e-yp(l + vo) 

K(r 2 , n) = - (p G h 2 - Ip ■ n 2 ) (3) 

where p = fq — ri, p = |p|, p = pjp, ® denotes the dyadic 
product, h 2 is the outward normal to the surface at fq, 
and I is the unit tensor. The expression for the Casimir 
energy is 


£ = — {^{iy) - ^{ioo)) dy (4) 

Jo 

where 

Hw) = '^rniiy) (5) 

even m'^2 

The functions are given by 

= q^-^TrK(r„,r,„_i)...K(f2,ri)K(ri,r^) (6) 

2 m dy 

where the operation Tr includes the integration of the 
variables ry over all surfaces, as well as the trace of the 
product of the tensors K(ri,ri_i). The terms can be 
represented as closed paths that visit m sites on the sur¬ 
faces. The paths having an odd number of points make 
no contribution because reversing the path changes the 
sign of each factor K, and both paths occur within the 
integral over paths. 

The sum in Eq.([^ can be done, giving 

q>{iy) = -y-^Tr ln(l - K • K) (7) 

dy 

This is the starting point for many calculations of the 
Casimir interaction energy ElIZ]. 

III. CASIMIR INTERACTION 

We have shown elsewhere [16] how to use this result 
to calculate the self-energy of a single conducting object. 
It applies equally well to the calculation of the energy of 
interaction of two objects: each surface integral is over 
the union of the surfaces of the two objects, and can be 
expanded into the 2"* assignments of the points of the 
paths to one surface or another. 

Several important simplifications result. The first is 
that we can drop the paths that are entirely confined 
to one surface or another, because these give the self¬ 
interaction of one of the two objects and is not an in¬ 
teraction between them. Although there are relatively 
few of these paths, they are the most singular, because 
they can be arbitarily short; removing them guarantees 
a finite path length and thus removes the need for the 
subtraction —$(ioo) in E q.([4] ). This in turn allows an 
integration by parts in Eq.M, so that the expression for 



4 


the Casimir interaction reduces to Sint = Sm>i^ 2 m, 
with 

£2Tn = -he [ -^TrK(rm,rm_i)...K(r2,ri)K(fi,rm) 

Jo 27rm 

( 8 ) 

where the Tr again includes integrals of each of the 
over all surfaces, but excluding the cases where all points 
are on the same surface. 


A. Two-point term 


The final simplification is that for small m it is prac¬ 
ticable to do the y integration in Eq.Q first. For the 
two-point contribution the only relevant paths go from 
surface of conductor A to that of conductor B and back 
again, with the result for the two-point function 





{p-dSA){p-dSB) d fe-^yp\ 
P dp\p‘^J 


where p = f a —tb, r a [tb) is the point on surface A (B) 
and dSA (dSB) is the vector surface element at Sa (tb)- 
The interaction energy is then given by 


S 2 = 


5hc f {p ■ dSA){p ■ dSB) 


IGtt^ 


P' 


( 10 ) 


The interaction (10) depends on the relative orientation 


of the surface elements; specifically anti-parallel surface 
elements attract while parallel elements repel each other. 
Application of Gauss’s theorem to Eq.(lO) once and 


then once again gives two further representations of the 
result 


see, Eq.(lO) is a good approximation for many standard 


problems. 

If the distance D between the two conductors signif¬ 
icantly exceeds their linear sizes, the interaction energy 
(121 can be written explicitly as 


£■2 = - 


IbhcVAVB 


(13) 


which agrees with the functional form of the Casimir- 
Polder original result m- For two spherical conductors 
of radii A and B we find £2 = —20hcA^B^/{3TTD'^). 
The exact result in this case m HO], Sint — 
— lA3hcA^B^/{16TTD'^) has about 25% larger amplitude. 

Similarly, for interaction energy of conductor A placed 
a large distance D away from conductive plane, Eq.(ll) 
predicts 


p _ 3hcvA 


(14) 


This again agrees with the functional form of the 
Casimir-Polder result m- For spherical conductor of 
radius A we find £2 = —hcA^/2TTD^. The exact result in 
this case [50], Sint = —QhcA^/lQ-KD^^ is larger in magni¬ 
tude by a factor of 9/8. 

According to Eq.([I^, the two-point approximation to 
the Casimir energy can be interpreted in terms of a pair¬ 
wise scalar interaction between uniform density clouds of 
’’Casimir charge” that fill the two objects. By doing the 
integral over one of the objects we can define a ’’Casimir 
potential” VB(r) due to conductor B such that the two- 
point contribution to the energy is given by the integral 
of the potential over the interior of the second object A-. 


£2 = j VBir^dvA 


(15) 


£2 = - 


15hc f {p ■ dSB)dvA 
IGtt^ / 


£2 = - 


15hc f dvAdvB 


47r3 


P' 


( 11 ) 


( 12 ) 


where dvA (dvB) is the volume element for conductor A 
(B) and the corresponding integrals are over the interior 
of one conductor or the other. Eq.(12) may be taken as 


an indicator of approximate additivity of Casimir-Polder 
interactions [H]. This was assumed by Mostepanenko 
and Sokolov m who proposed the functional form 
to complement the PFA. We see that the BD method 
justifies the conjecture m and derives an amplitude; in¬ 
cluding higher order terms allows for systematic improve¬ 
ment on Eqs.( l0p2 ). 

Eq.(12) implies that interaction is necessarily attrac¬ 


tive (to this order of approximation) for any pair of ob¬ 
jects, so that a repulsive Casimir interaction will require 
that the higher order terms play a dominant role. Eqs.( 10 


12) are equivalent representations of the first term of 


an expansion which we believe to be exact. As we shall 


For simple geometries this can be used to find closed 
form of the interaction energy at arbitrary separation of 
the conductors. For the plane, the potential depends only 
on the distance z from it, in the form 


Pp/ane(2') — 


She 


(16) 


• Integrating this over the half-space z > D gives the 
interaction per unit area of parallel planes. To this order 
of approximation 


£2 he 

A 


(17) 


which differs from Casimir’s result m Sint j A — 
-7r2;ic/(720D3) by a factor C(4) = tt^/QO = 1.0824 (BD 
[ 5 ] show how the higher m contributions £m give rise to 
this correction). 


Integrating the potential (16) over a sphere of ra¬ 


dius A at distance D from the plane gives the interaction 
energy for sphere and plane 


£2 = - 


heA^ 


27r(Z)2 - A2)2 


(18) 
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For D ^ A this agrees with Eq.(14) while for D 
becomes [5] 


£ 2 -- 


hcA 


he 


hcA 


-I 


A it 


(19) 


87rX2 SttX 87rX2 
where X = D — A \s the distance of closest approach. 


The first term in (19) is very close to the PFA result 


(again too small by the factor C(4)) while the second 
represents a correction to the PFA in terms of the di¬ 
mensionless proximity parameter XjA. It compares rea¬ 
sonably well with the result of an asymptotically ex¬ 
act calculation of the sphere-plane interaction energy 
£mt{X ^ A) oc (1 — 1.69X/A) [T3] especially given that 
Eq.([T^ is just the two-point result. Below we will show 
that it also works reasonably well at all separations. The 
same functional form as in Eq.(18) (with a different am¬ 
plitude of phenomenological origin) was given previously 

m- 

For a sphere of radius A, the potential depends on the 
distance r from its center 


^sphereij'^ — 


hcA^ 5r^ -I- A^ 

TT^r (r2 — A2)4 


( 20 ) 


For r « A this reduces to the potential of the plane (161. 


The interaction energy with a second sphere of radius B 
and center at D is 


£2 = - 


hc 



R 12 R 21 

R 11 R 22 

AB 

AB 

Wi 

+ ^ 
^12 


2D2 


2D2 


R 11 R 22 R 12 R 21 


+ 


AB AB 


Rh 


p2 

-^22 


( 21 ) 


where Rn = D A 
D — A + B, and R 22 


T B, Ri2 — D 
= D - A - B. 


-j- A — B, R 21 — 
For large D this 


reduces to Eq.(13). For D close to A -|- i?, we find £2 ~ 
-hcAB/[8Tr{A + B){D -A- B)^]. 

Other geometries that do not allow explicit analytic 
calculation of £2 can be studied by numerical evaluation 
of Eqs.([To|), ( [IT] ), or (fl^. We observe that when the 


separations between objects is large, Eq.(10l is less ap¬ 


propriate because the near and far sides of each object 
have contributions of opposite sign; when the objects are 
near each other, however, Eq.(lO) may be more efficient. 


• These encouraging results are offset by geometries in¬ 
volving macroscopically long cylindrical conductor. For 
example, for a cylinder of length £ and radius A at dis¬ 
tance D ^ A from conductive plane, we find for an in¬ 
teraction energy per unit length 


£2 ihcA^D 

Z “ ~ 87r(£)2 - A2)5/2 


( 22 ) 


The same functional form (with different amplitude) was 
given previously m- For D A this agrees with 
Eq.(14| while for « A it becomes 


3feAi/" 1^1 

^ 32y27rA:5/2 4 ^ a ) 


( 23 ) 


where the leading order term is very close to the PFA re¬ 
sult [21] {£ 2 /£pFA = C~^(4) = 90/Z = 0.9239) while the 
PFA plus proximity correction, £2 oc (1 — XjAA), com¬ 
pares reasonably well with asymptotically exact result, 
£int oc (1 — 0.48103A1/A) [T^]. However, at intermediate 
cylinder-plane separations D :s> A (where £ is the largest 
length scale in the problem), both Eq.(|2^ and the PFA 
are qualitatively different from the result due to Emig et 

al. ED, £^nt/£ ^ -hc/{D^ HD/A)). 

Similarly, for the Casimir interaction of two well- 
separated cylinders of arbitrary cross-sectional areas a a 
and (Jb find 

Z _ dfecTAO-B 

£ 7r3£»6 ^ 

which must be valid if D is the largest length scale of 
the problem. However, this disagrees with the asymptot¬ 
ically exact result due to Rahi et al. E3 cylinders 
of radii A and B, £intl£ — —hc/[D^hY{D/A)hY{D/B)\ 
valid when D ^ A, B and £ is the largest length scale 
of the problem. These anomalous cases will be further 
discussed below. 


B. Four-point interaction 


The dependence of K, Eq.([^, on the imaginary 
wavevector y is entirely contained in the scalar prefac¬ 
tor, so that the integral over y in Eq.(|^ can always be 
done first. For the four-point term £4 the integrand is 

pOO 

•^= (1 + 2/P43)(1 + yP32)(l + 2 /P 21 ) 

Jo 

X {l + ypi 4 )exp{-y[p])dy (25) 


where the abbreviation [p"] = pjg -|- P 32 -I- P 21 + P 14 fi^-® 
been introduced. A little algebra shows that 


J = 3[pH-[pHH] 


6 [p][l/p] + 24 


(26) 


in which we have introduced the further abbreviation 
n = P 43 P 32 P 21 P 14 - This reduces the problem to the eval¬ 
uation of 


Jic 7 

£4 = (2^)ZP 2)k(2, l)k(l, 4) (27) 


where k(2,1) = p 2 i ® H — Ip 2 i • H is the tensor part 
of K, and the Tr includes an integration of each r( over 
both surfaces, excluding the cases where all four points 
are on the same surface. For all but the case of parallel 
planes, this expands out into a large number of difficult 
integrals, and analytic progress seems unlikely. However, 
the numerical integration is straightforward. 


C. Higher-order contributions 

This process could be continued to higher m. The an¬ 
alytic integration over y can be done as described above. 
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but the result will be algebraically complex. Since we 
will be forced to use numerical integration over the 2 m 
two-dimensional surfaces (an integration over a 4 to di¬ 
mensional space), we can just as well do the numerical 
integration over y at the same time. We will show be¬ 
low that for selected problems, the first two £m give us¬ 
ably accurate values. We will also discuss cases involving 
cylinders, where carrying on to small finite order is in¬ 
sufficient. 


D. Comparison with other approaches 

Most routes to the evaluation of the Casimir interac¬ 
tion have made use of a contour integral representation. 
The integrand of the contour integral contains a response 
function that has a pole at every mode frequency. In the 
present case, this response function is the relationship 
between js and j implicitly defined by Eq.([^. 

Rather than directly solving ([^, BD represent the so¬ 
lution as an infinite series by repeatedly substituting the 
equation into itself. This has the advantage that the ex¬ 
pression for the Casimir interaction is in the form of a 
sum of integrals over known functions. 

Analytic progress can be made on the direct solution of 
([^ in cases of high symmetry, where the tensor K has a 
diagonal representation. This has made possible the eval¬ 
uations for the interaction of two spheres mus], sphere 
and plane [23ll25] . cylinder and plane [21], and parallel 
cylinders [H]. However, there is little further that can be 
done by this approach, since even for the prolate ellipsoid 
it is unknown how to separate Maxwell’s equations |26j . 
Each of these evaluations requires its own mathematical 
elaboration (expansions in Bessel functions or spherical 
harmonics), which allow very little generalization. 

Subdividing the surfaces of interacting objects into 
small polygons converts into a large matrix equa¬ 
tion, which can be inverted by standard algorithms. This 
approach has been explored by Johnson and collabora¬ 
tors |8] . Although this does afford a general route to the 
evaluation of the interaction between objects of arbitrary 
shape, its use requires some expertise in the art of tes¬ 
sellation of surfaces and inversion of large matrices. The 
output of the program is the energy of one configuration 
(or the stress one exerts on the other); getting a sense for 
how the objects interact can only come from studying the 
numerical output. 

In contrast, the BD expansion describes the interac¬ 
tion in terms of an integral ([^; we have shown that the 
integral over imaginary frequency can be done, leading, 
in the two-point approximation, to the representations 
0, (HI, whose form is simple enough that ana¬ 
lytic evaluation is possible for the kinds of problems that 
have been worked in the past, and that one can learn 
something about the nature of the interaction from it 
even when analytic evaluation is not possible. The two- 
point approximation is a controlled improvement on the 
PEA: about as easy to use and much more general. 


We have shown that the integral over imaginary fre¬ 
quency can also be done for the four-point interaction, 
leading to (27), which now has the form of a geometrical 
interaction between the objects. As will be shown below, 
the four-point interaction often is a small correction (in¬ 
dicating rapid convergence of the BD series). We discuss 
the exceptions below. 

In the forms that we have presented them, the two- and 
four-point interactions can be readily evaluated numeri¬ 
cally. We have chosen to use a rather basic implementa¬ 
tion of the Monte Carlo method (averaging the integrand 
over randomly chosen multiplets of points) because it is 
easy to code, and the same set of programs could be used 
for all of the geometries we describe, replacing only the 
subroutine that describes the interacting surfaces. Once 
it had been set up, the program took about a day (run¬ 
ning on a standard desktop computer) to generate the 
figures we display. 


IV. APPLICATIONS 

To judge the utility of the BD expansion, and in par¬ 
ticular the accuracy of the two-point term ( |Io| regarded 
as a replacement for previous approximations pun] , we 
consider some examples. 


A. Parallel planes 

Problems in which one of the surfaces is a plane are 
simplified because the trace of the product of tensors K 
appearing in (|^ vanishes when two successive points are 
on the same plane. Eor parallel planes this implies that 
the integration points must alternately be on one surface 
or the other. BD discuss this problem and show that the 
successive contributions differ only in a factor of 
explaining why the two-point contribution is less than 
the exact answer by the factor C(4) = = 1.0824. 


B. Plane or sphere facing arbitrary surface 


Eq. (10) simplifies quite a bit when surface i? is a plane 


facing surface A whose height relative to B is given by a 
function z(x,y). The integration over the plane can be 
done to give 


Thus for this problem, the first term of the BD expan¬ 
sion (up to the factor of C(4) = tt'^/OO = 1.0824) derives 
the PEA; further corrections come from the higher order 
terms. Our earlier results for the interaction of the plane 
with the sphere, (18), or the cylinder, (22), can be repro¬ 


duced as a superposition of the interactions with front 


and back surfaces computed with the help of Eq.(28). 
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Figure 1: (Color online) Casimir interaction energy for a 
sphere of radius A with center at distance D from a plane 
(sketched in the inset). Green line: Results from Ref. |23| . 
Red line: Results from Ref. |24| . Blue line: sum of two- and 
four-point approximations. These have been divided by the 
two-point result, Eq.( |18| l, which removes most of the depen¬ 
dence on D (the straight black line represents the two-point 
approximation itself). At the right we have indicated the 
short distance asymptotic value £int/£-i = 7r^/90 = 1.0824 
due to Casimir [1] as discussed in Section IVB. At the left we 
indicate the large distance asymptotic value Eint/S -2 = 9/8 = 
1.125 due to Boyer |20 |. 


A similar reduction of the integral over both surfaces 
to an integral over only one can be effected for the case 
of a sphere of radius A and an arbitrary surface R{9,<j)) 
in the form 


£2 


Ahc f 5R^-A^ 

“ 4 ^ J °^‘^^(i?2-A2)3 


(29) 


C. Sphere and plane 


Emig [24], Kenneth and Klich |23| and Maia Neto, 
Lambrecht and Reynaud |25j have evaluated the inter¬ 
action energy for a sphere of radius A with center at 
distance D from a plane. 

Figure 1 compares the results of Emig and those of 
Kenneth and Klich to the sum of the two- and four-point 
terms, in each case dividing by the two-point result ( [I^ . 
The Casimir interaction varies by six orders of magnitude 
over the range of the data in this Figure, but the two- 
point approximation captures all but 10% of this. The 
four-point term accounts for most of the discrepancy. 


Figure 2: (Color online) Casimir interaction energy for 
spheres of radius A at distance D from each other (sketched 
in the inset), divided by the two-point interaction energy £ 2 - 
Black straight line: The two-point interaction Eq.(30L Blue 
line: includes the four-point interaction. Red line: Results 
from Ref. [23] . At the right we have indicated the asymp¬ 
totic value Sint /£2 = 7r'^/90 = 1.0824. At the left we indicate 
the asymptotic value £intl £2 = 429/320 = 1.34. 


D. Interaction between spheres 


The interaction between identical spheres of radius A 
and with centers separated by distance D has been stud¬ 
ied by Emig et al. [7] and by Kenneth and Klich [23j . 
The two-point approximation is the special case A = B 


of ( 21 ): 


^ _ lie / / £>2 \ 40A4-6A2D2 2^2] 

+7S2'%^4';42)2 " + J 
(30) 

The four-point interaction was evaluated by constructing 
four points Ty randomly chosen to be somewhere on ei¬ 
ther surface. The cases where all four were assigned to 
the same surface were discarded; for the rest, the inte¬ 
grand for £ 4 , Eq.(|^, was evaluated and averaged over 
10® attempts for each value of D. 

Figure 2 compares £2 + £4 to the results of Ref. [23] 
(which are nearly identical to those of Ref.jTj). Since 
these are strong functions of D, we have divided each 
result by £ 2 - 


E. Hyperboloid and plane 

The surface z{x,y) = D + A A — o?A 

describes a cone with a rounded tip, as might be used 

























1.2 


SFS 


SBS 


SFBFS 



Figure 3: (Color online). Casimir interaction energy between 
a hyperboloid and a plane (sketched in the inset), as a func¬ 
tion of distance of closest approach D, for various asymp¬ 
totic opening angles. Blue: a = 1.5; Red: a = 1.0; Green: 
a = 0.75. 


in an atomic force microscope. The parameter a is the 
cotangent of the opening angle of the cone; the length 
scale A is the radius of curvature of the tip; and D is 
the distance of closest approach to the plane z = 0. Eq. 
(281 can be used to find the interaction in the two-point 
approximation 


£2 = - 


hc{D -|- a^A) 
Sna'^D'^ 


(31) 


For D ^ a^A this reduces to the small-distance limit 
of the interaction between sphere and plane, as given 
by Eq.(18) (after making the replacement D ^ D + A 
in that equation). The limit a —>■ 0 cannot be taken, 
because the total interaction energy for parallel planes is 
infinite. In the a —>■ 00 limit the hyperboloid turns into a 
paraboloid z{x,y) = D + {x‘^ +y‘^)/2A whose interaction 
energy with the plane in the two-point approximation is 
£2 = -HcA/{8TrD‘^). 

We evaluated the four-point term as described above, 
with the results shown in Figure 3. As previously, we 
have divided the sum of £2 + £4 by £ 2 , Eq.(31), to facili¬ 
tate comparison. 


F. Sphere and slab: illustration of screening 


We can appreciate the need for higher-order terms in 
the BD expansion by considering the case of a sphere 
of radius A at distance D from an infinite plane slab of 
thickness H. According to (101 the front surface of the 
slab will attract the sphere, while the back surface repels 
it; according to (12) the net interaction is attractive but 



Figure 4: (Color online). Top: Conducting sphere facing a fi¬ 
nite thickness slab (shown in grey); two-point (SFS and SBS) 
and four-point (SFBFS) scattering diagrams are sketched. 
Bottom: contribution of these diagrams into the Casimir in¬ 
teraction energy; the four-point diagram SFBFS essentially 
cancels the effect of the back surface represented by the SBS 
diagram thus illustrating screening. Numerical data employed 
to produce the graphs are assembled in the Table. 


will become small for small H. In the limit D > A, 

= -^ (i^ - (idw) 

However, the interaction of the sphere with a conducting 
plane is attractive, with no need for consideration of what 
might lie beyond it. This contradiction arises because the 
two-point interaction describes an unscreened interaction 
between the surfaces, yet in the envisioned situation the 
front surface of the slab will greatly screen the interaction 
between the sphere and the back surface of the slab. This 
screening is described by the higher order terms of the 
BD expansion. 

The two terms of ( [3^ arise in (1^ by the cases where £2 
is on the the front (Fjor back (B) surfaces, respectively, 
while rq is on the sphere (S). Drawing these as diagrams, 
sketched in Figure 4 (top) the first term is a closed path 
from sphere to the front surface of the slab and back 
to the original point on the sphere (represented as SFS). 
The unwanted interaction is a closed path from sphere to 
the back surface and back to the sphere (SBS), crossing 
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the front surface both ways. The four-point interaction is 
represented by a large number of diagrams, most of which 
are corrections to the terms already calculated (such as 
SSFSS) or identically zero. We believe that the relevant 
case is SFBFS, where the interaction between sphere and 
back surface is modified by crossing the front surface. We 
calculated this for a unit-radius sphere and unit-thickness 
slab (^ = 1, i? = 1) by a Monte Carlo integration as 
described in Section HID. To simplify comparison, we 
have divided all results by SFS = —CjD'^ (a proxy for 
the right answer), with the following results: 


D 

1 -b SBS/SFS 

SFBFS/SFS 

sum 

error 

SBFBS/SFS 

5.000 

0.530 

0.468 

0.998 

0.005 

-0.274 

10.000 

0.319 

0.696 

1.015 

0.028 

-0.501 

12.000 

0.275 

0.786 

1.062 

0.041 

-0.562 

15.000 

0.228 

0.838 

1.067 

0.129 

-0.661 

18.000 

0.195 

0.783 

0.978 

0.141 

-0.622 

20.000 

0.178 

0.681 

0.859 

0.183 

-0.925 


The fourth column shows that SFBFS has substan¬ 
tially canceled SBS for all D (the fifth column is an esti¬ 
mate of the statistical error of the Monte Carlo integra¬ 
tion). We note, however, that another of the four-point 
diagrams, SBFBS, is large. We believe it will be canceled 
at higher order by diagrams such as SFBFBFS. The data 
shown in the Table are also displayed in Figure 4. 

The considerations of this paragraph would also be rel¬ 
evant to other systems involving large, thin objects, such 
as an extremely oblate ellipsoid with a sphere near its 
axis. They also explain the success of the PFA (which 
ignores the effect of the back surface) at short separation 
distances. 


G. Configurations involving cylinders 

We now assume that in the plane-cylinder geometry 
the cylinder length H is effectively infinite, i.e. the inter¬ 
action energy is strictly extensive in the cylinder length, 
and focus on the large separation regime D ^ A w her e 
the two-point prediction, E 2 IC ~ —hcA'^/D'^, Eq.([^, 
is in qualitative contradiction with the exact result j21j . 
E,r.tlC^-hcl{D‘^ HD/a)). 

1. Semi-quantitative argument 


A similar argument applied to the problem of a spherical 
conductor of radius A placed a distance D away from 
conductive plane predicts 


Ein.t — 



(34) 


Here fc,sit) are dimensionless functions determined by 
the geometry of the problem. Both in the case of the 
cylinder and sphere one must have /c,s(0) = 0 because 
absence of the conductor implies absence of interaction. 
However, the character of the approach fc,sH/D —>• 
0) —0 is qualitatively different in the two cases. Intro¬ 
duction of an infinitesimally small spherical conductor 
is a small perturbation. Therefore its interaction with 
the plane as A/D ^ 0 must be proportional to the con¬ 
ductor’s polarizability, i.e. to its volume, thus implying 
fs{t —>■ 0) (X This leads to the Casimir-Polder type 
result ( |l4| . 


Introduction into space of an infinitesimally narrow 
cylindrical conductor, on the other hand, is a singu¬ 
lar perturbation because distribution of the conductive 
material in space is not compact even if A is infinites¬ 
imally small. The difference between the sphere-plane 
and cylinder-plane geometries is topological - in the for¬ 
mer case any contour in vacuum by continuous defor¬ 
mation can be contracted into a point while this is not 
true in the latter case. For the cylinder-plane geometry 
the topology of space changes at A = 0. Therefore the 
dependence fc{A/D —)• 0) must be non-analytic. 

If topology is indeed the fundamental reason behind 
the anomalous dependence on separation distance, then 
non-analytic dependence of the interaction energy on the 
vicinity to the point of topological change will be rele¬ 
vant to any problem for which the distance between the 
attracting objects is not the largest length scale. 

In order to determine the dependence fc{A/D) we first 
look at a simpler problem which belongs to the same 
universality class: the cavity formed between the two 
concentric infinitely long cylindrical conductors p > D 
and p < A, i.e. a coaxial cable. The reference state here 
will be chosen to be the absence of the inner conductor, 
A = 0. Then the interaction energy due to insertion of 
the inner cylinder of radius A, directly given by Eq.Q, 
must have the functional form (33). 

The TEM mode should be omitted as it does not have 
limiting frequency set by spatial scales of the problem. 
The dispersion law of the TE and TM modes is given by 


Anomalous interactions arising when long cylinders are 
involved can be understood via a combination of dimen¬ 
sional analysis and physical reasoning. In the limit that 
the length C is large compared to other dimensions, we 
discuss the interaction energy per unit length, which de¬ 
pends on the parameters h and c, and on two length scales 
D and A. Dimensional analysis implies that 



(33) 


H = c^{q^ -b K^) (35) 

where q is the continuous wave number along the cylin¬ 
der axis while k represents the discrete set of transverse 
wave numbers corresponding to the cutoff frequencies. 
The summation over a in Eq.([^ now stands for a sum 
over TM and TE modes with integration over q and sum¬ 
mation over all possible k. 

The spectrum of the TM modes is given by solutions 
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to the transcendental equation 


_ Jn{nA)Nn{KD) 

Nni^A) 


(36) 


where J„ and Nn are the Bessel functions. The spectrum 
of the TE modes is given by similar equation except that 
all the Bessel functions are replaced by their derivatives. 

At ^ = 0 both the TM and TE equations reduce to 
JninD) = 0 and J'^{kD) = 0 which give us the spec¬ 
trum of the cylindrical cavity of radius D. We are now 
going to follow solutions to these equations as small A 
is introduced. We already know that A = 0 is the point 
of non-analyticity of the energy, so we need to look at 
displacements of the transverse wave numbers which are 
non-analytic in A. Employing the small argument expan¬ 
sion of the Bessel functions we then infer that only the 
n = 0 TM mode produces shifts that are non-analytic in 
A. So we need to focus on the equation 


ri = {Acos (pi, AsiTKpi, Zi) are on the cylinder, making 
approximations appropriate to the limit D ^ A. These 
paths that make only one visit to the plane are sufficient 
to explain the result, and are the most important contri¬ 
butions to ([^ owing to the exponential dependence on 
path length. Note that in the surface integrals in Q each 
of the points could be on the plane, and that relabeling 
so that the point on the plane is always r 2 eliminates the 
factor of 1/m in Q. Also note that ([^ vanishes when 
both fq and fq are on the plane. 

In cylindrical coordinates, a current zjc{Z) (indepen¬ 
dent of 4 >) on the cylinder gives rise to a magnetic field 


BoutipAjZ) = (p dZ' 


7 °° ^^'‘‘^^~^'^kAIo{kA)K'o{kp)j,{Z') 

(40) 


Jo{x) = 


Jo(xt)No{x) 

No{xt) 


X = kD, t = 


A 

D 


(37) 


This equation is solved iteratively around solutions to 
the A = 0 problem Jq{xi) = 0 where I labels the roots 
of the zero-order Bessel function. The lowest order in 


A/D << 1 solution to (37) for the displacements Axi 
due to hnite A is 


Ax, 


ttNo(xi) 


2 J^{xi)HA/D) 


(38) 


Knowing all the displacements of the transverse wave 
numbers Ani = Axi/D one can find the change in the 
energy produced by having A small but finite. However 
as far as the A-dependence is concerned, this is where 
we can stop because the entire dependence on A is given 
by the multiplicative factor of 1/ ln(A/i7), thus implying 
fc(t) oc 1/lnt. The remaining dependence of the interac¬ 
tion energy per unit length on D follows from dimensional 
analysis, Eq.(33): 


£i' 


he 


c 


£>2 ln(£)/A) 


(39) 


We see that anomalous dependence on the separation 
distance is due to the n = 0 TM mode. 


2. Summation of the multiple seattering series 

We will now show that for the cylinder-plane geometry 
the same anomalous dependence on separation distance 
(39) follows from a summation of the entire BD series, 
and that it is a generic feature of the Casimir effect for 
geometries involving cylinders or wires. 

The geometry is a cylinder of radius A along the 
Z axis and an infinite conducting plane a distance D 
from it along Y. We will evaluate for the case 
that fq = {X,D,Z 2 ) is on the plane and all other 


for p > A and by 


Bin{p,4>,Z) = p dz' 

^ —<50 

X J ^e*«(^-^')fcAKo(fcA)/'(fcp)jc(^') 

(41) 

for p < A, where Iq and Kq are the modified Bessel 
functions and k'^ = +y‘^. The field is discontinuous at 

the cylinder surface because it is carrying a current; there 
also is a field due to currents elsewhere (on the cylinder 
and on the plane). This nonlocal part of the field is 
canceled by the local field inside the cylinder, according 
to the rule 


jc{Z) — 2p X -{Bout + Bin) 


(42) 


Substituting (40) and gives a relationship between 
the current at different points on the cylinder in the form 
of an integral similar to ([^ involving a kernel Kcc(ri, 7*2 )■ 
We can take advantage of the translational symmetry 
along the cylinder by using a Fourier representation. De¬ 
fine 


According to the argument above, the explicit form of 

Woe is 

W,, = fcA[/o(fcA)K'(fcA) + /'(fcA)Ko(fcA)] 

= l-2AkIo{kA)KoikA) 

~ 1- k^A^lnil/kA) (44) 

where the small argument approximations for the Bessel 
functions are taken at the last step; those are appropriate, 
since we will only make use of this result for kD < 1 and 
A D. By the same argument we will replace ln(l/fcA) 
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by \x\{D/A). This argument can be repeated for the case 
that the current on the cylinder varies as cos n(j). The 
effect is to replace the Bessel functions of order 0 by 
functions of order n. The important difference is that in 
the small k limit the corresponding Wcc is well less than 
unity. In the argument to follow, this will imply that the 
cylinder does not have resonant response to this part of 
the current, so that it can be neglected. 

The field (401 induces a current j{Z,X) = 2Y xB, and 
this current gives rise to a field that can be calculated 
similarly to (41) and which in turn requires a screening 


field on the cylinder. The corresponding kernels K can 
also be given in a Fourier representation, although these 
also depend on the coordinate X of the plane: 


K(fi,r2) = (45) 


The explicit forms of the response functions are 

2ADk 


Wp,(<7,X) = - 


R 


-Io{kA)Ko'{kR) 


W,piq,X) = ’^I'^{kA)Ko{kR) 

7T 


(46) 


where R = \/X'^ + 

After substituting this into , all of the surface inte¬ 
grals except the integration over X can be done immedi¬ 
ately, leading to 



X 



dq 

27r 



Wep(X)Wp,(X)dX ^ [Wee]'"-' 
m>2 




dX 


dy r ^ 

2 tt 2tt 

Wep(X)Wpe(X) 
1- Wee 


(47) 


Now the significance of the small argument behavior of 
Wee becomes clear: because it approaches unity, an ar¬ 
bitrarily large number of terms of the expansion must be 
included. The cylinder is exhibiting a resonant response 
to the presence of the plane. 

Substituting the explicit forms for the response func¬ 
tions, the integrals can be done to reproduce the result 

of Ref. [21]: 


~167rD2 ln(i:»/A) 


(48) 


valid for H ^ A. 

From this discussion it is clear that a similar anomalous 
Casimir interaction will occur for any geometry involving 
a cylinder or wire whose length significantly exceeds its 
radius of curvature. 


V. EFFECTS OF FINITE TEMPERATURE 


Ideas described up to this point can be extended to 
the case of finite temperature when the formula for the 
Casimir energy is replaced by an expression for the 
Casimir free energy |^: 


he 

T = — / ($(zj/) — $(foo)) dy 

TT do 

+ 2T r^^[$(*y)-$(+*0)]y(y) 
do y 


(49) 


The temperature effects are accumulated in the second 
integral and the temperature T is also present in the 
sawtooth function 


, . 1 he 

^^y^ = 2 -^y 


Eefd- 


2'nnT \ 

he J 


(50) 


where 0(a:) is the step function. Similar to the zero- 
temperature case, Eqs.([^ and can be substituted into 
Eq.(|4^ so that the expression for the Casimir free energy 
of interaction reduces to Xmt = J^m>i with J' 2 m 
standing for the 2TO-point contribution. Since at T = 0 
and many geometries the two-point contribution captures 
the bulk of the effect, we expect the same to hold at finite 
temperature; some of the evidence to that is given below. 
In what follows we limit ourselves to the two-point con¬ 
tribution T 2 into the Casimir free energy of interaction 
of conductors A and B. Combining Eqs.([^,(49l and (50) 
we find 


T 2 = 



{p-dSA){p-dSB) 


1 , 2'KpT \ 

- coth —^ 
p Tie J 


(51) 


Compared to the zero-temperature result (10), the novel 


feature of Eq.(51) is the presence of cross-over scale 


A = 


he 

2 ^' 


(52) 


the well-known Casimir length, that defines two regimes 
of separation of the conductors: 

(i) For p <C A thermal fluctuations are negligible and 
the zero-point motion dominates the physics. This is the 
regime of Casimir/Casimir-Polder interactions. At T = 0 
when A = 00 Eq.(51) reduces to the zero-temperature 
result 


( |10| ). At T = SOOiF we find A = 1pm, thus imply¬ 
ing that even at room temperature quantum mechanics 
dominates Casimir interaction of sub micron separated 
conductors. 

(ii) For p 3> A thermal fluctuations dominate the 
physics and zero-point motion is negligible. This is the 
regime of Van der Waals interactions. In the high- 
temperature regime the Planck’s constant drops out of 
Eq.(51) which acquires purely classical form 


^2 = 


r {p-dSA)ip-dSB) 

47r2 / p6 


(53) 
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As in the zero-temperature case, this can be interpreted 
as a superposition of vector interactions of elementary 
surface elements; in contrast to the T = 0 situation the 
interaction (53) is more long-ranged. Application of the 


Gauss’s theorem once and then once again gives two fur¬ 
ther representations of the result 


T f{p-dSB)dvA 


^2 =- 


3T f dvAdvB 
27r2 / n® 


(54) 


(55) 


which are high-temperature counterparts of Eqs. 0 and 
(12), respectively. Eqs.( 53p5 ) have the form J 2 = —TS 
where S is the Casimir interaction entropy thus imply¬ 
ing that in the high-temperature regime the Casimir ef¬ 
fect has purely entropic origin. Eq. (55) also tells us that 
5 > 0, i.e. the interaction is attractive (to this order of 
approximation) for any pair of conductors. 

Eor simple geometries (plane-plane, plane-sphere, 
sphere-sphere, and plane-cylinder) it is possible to derive 
closed form high-temperature counterparts of the results 
of Section IIIA. Eor example, for two conductors sepa¬ 
rated by a large distance D we find 


^2 = - 


3Tvavb 

27r2H6 


(56) 


which captures the functional dependence of the exact 
result |18| and has the same order of magnitude. 

Similar to the analysis of the zero-temperature and 
high-temperature limits the general expression (51) can 


be transformed via Gauss’s theorem to find somewhat 
cumbersome formulas encompassing Eqs. ([g, (@, 
and (55) as special cases. The main physics consequence 


is that interaction of any two conductors is attractive (to 
this order of approximation) for arbitrary temperature. 
Eor the two-plane geometry the integral (51) can be com¬ 


puted for arbitrary temperature, and the interaction free 
energy can be given in closed form. This would repro¬ 
duce the very accurate result due to BD [S] who found it 
differently by specializing the two-point expression (|^ to 
the case of the two-plane geometry. Closed form expres¬ 
sions for the interaction free energy can be also given in 
the Casimir-Polder limit, i.e. for the interaction between 
well-separated small particles and between small particle 
and the plane. The outcome, which is generalization of 
Eqs. (13) and (14) to the case of finite temperature, has 


the same functional form as the result of complementary 
calculation m given in terms of particle polarizabilities. 
Other geometries that do not allow explicit analytic cal¬ 
culation of T 2 can be studied by numerical evaluation of 

Eq.(|ni). 
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